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onboard  targeting  computations. 
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ABSTRACT 


Methods  are  investigated  which  apply  to  the  conversion  of  latitude,  longitude 
and  height  to  cartesian  coordinates  in  a  plane  tangent  to  the  earth's  surface  in 
connection  with  onboard  missile  targeting.  The  criterion  for  the  usefulness  of 
the  method  is  the  error  in  the  north  and  east  coordinate  with  reference  to 
Clarke's  spheroid  1866.  This  error  is  determined  as  a  function  of  vet. tor  length 
and  a/imuth  for  the  spherical  earth  model  referenced  to  the  mean  geodetic 
latitude  at  White  Sands  Missile  Range  utilizing  ?.  program  prepared  for  a  pro¬ 
grammable  electronic  desk  calculator  (Marchant  1016  PR).  An  approximation 
of  geocentric  latitude  is  used  in  the  program.  It  is  explained  ho’v  the  error  in¬ 
herent  to  the  spherical  earth  model  car.  be  reduced  applying  a  certain  correction. 
An  expression  for  computing  the  reference  coordinates  in  a  plane  tangent,  to  the 
earth  spheroid  is  derived  which  requires  less  numerical  effort  than  the  standard 
p rocedure. 
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SECTION  I 


INTRODUCTION 

Approximations  for  Converting 
Geodetic  to  Cartesian  Coordinates 

The  purpose  of  a  study  related  to  the  SRAM  project  was  to  simplify  cal¬ 
culations  for  the  post-flight  evaluation  of  trajectory  data  and  to  reduce  th 
number  of  numerical  operations  so  that  a  programmable  desk  calculator  could 
be  used  for  the  computations.  Of  primary  interest  was  the  process  used  to 
convert  geodetic  da*a  of  latitude,  longitude  and  height  to  cartesian  coordinates 
in  a  plane  tangent  to  the  surface  of  the  earth.  The  Marchant  Calculator,  Model 
PR  1016,  was  available  for  these  calculations  and  it  was  required  that  position 
fixes  and  target  range  computations  performed  during  the  flight  of  a  B-52  air¬ 
craft  be  examined  on  the  ground  to  obtain  "quick-look"  information  revealing  the 
accuracy  of  the  onboard  computations.  This  calculator  had  to  be  utilized  to  save 

the  "turn-around"  time  spent  in  waiting  for  the  results  of  a  large  electronic 
computer  (CDC-3600). 

To  find  an  adequate  procedure,  conversion  methods  had  to  be  compared  with 
respect  to  their  accuracy  and  a  compromise  mude  between  the  complexity  of 
the  calculations  and  the  accuracy  of  the  final  results.  Known  approximations 
had  to  be  examined  and  their  limitations  with  respect  to  the  maximum  range  be  ¬ 
tween  the  fixed  point  or  target  and  the  aircraft  had  to  be  determined.  The  goal 
was  to  implement  the  calculator  with  a  program  capable  of  converting  geo¬ 
detic  data  to  tangent  plane  data  with  an  accuracy  equivalent  to  a  few  feet  com¬ 
pared  to  the  solution  of  using  a  spheroidal  earth  model  which  was  assumed  to 
yield  the  exact  solution.  A  specific  question  was  then,  up  to  what  maximum 
range  the  spherical  earth  model  would  be  adequate  in  this  respect. 

Target  ranges  in  the  White  Sands  test  area  barely  exceed  100  nautical 
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miles;  however,  because  of  the  significance  of  this  subject,  errors  caused  by 
some  approximations  were  calculated  for  ranges  up  to  700  nautical  miles.  These 
ranges  or  distances  between  the  aircraft  and  targets  are  calculated  to  determine 
the  "range -to-go",  which  is  derived  from  the  output  of  the  inertial  platform  sen¬ 
sors  and  the  known  target  coordinates. 

Since  information  about  approximations  for  the  required  conversion  process 
are  scattered  in  the  literature,  the  more  important  expressions  related  to  this 
subject  are  reviewed  and,  in  some  cases,  derived. 
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SECTION  II 


METHODS  FOR  THE  CONVERSION  PROCESS 

It  is  standard  practice  to  calculate  position  differences  of  an  aircraft  with 
respect  to  a  ground  target  in  cartesian  coordinates.  To  relate  the  two  positions 
to  each  other,  it  is  necessary  to  convert  the  geodetic  survey  data  of  the  target 
into  cartesian  coordinates  of  the  platform  system  or  vice  versa.  Normally  a 
north/east  orientation  is  used  as  a  reference,  and  a  deviation  of  the  platform 
axes  from  this  orientation  is  taken  into  account;  results  are  presented  with  com¬ 
ponents  in  north  and  east  direction  counted  positive  which  are  called  Rj^  or  X 
and  Rg  or  Y  hereafter.  Of  less  interest  here  is  the  Z  component  which  may  be 
computed  but  can  also  be  found  from  the  altitude  measurement  of  the  onboard 
radar  altimeter  combined  with  information  of  terrain  altitude. 

In  the  following,  it  is  assumed,  unless  otherwise  stated,  that  the  point  on 
the  earth's  surface  is  at  a  mean  sea  level.  The  models  investigated  with  respect 
to  their  accuracy  are: 

Circular  arcs  in  north/ south  and  east/west  direction  and 

Spherical  earth  with  geocentric  radius  derived  from  the  mean 
latitude  of  the  two  points. 

The  spheroidal  earth  model  with  major  and  minor  axes  specified  for  the 
basic  ellipse  serves  as  a  reference  to  determine  the  error  of  the  results  de¬ 
rived  both  from  circular  arcs  and  spherical  earth.  It  is  noted  that  the  spheroidal 
earth  model  is  also  an  approximation,  but  it  best  approximates  the  actual  earth 
shape. 

a.  Circular-arc  Approximations 

Circular  arcs  can  be  approximated  using  either  geodetic  or  geocentric 
quantities  for  the  earth  radius  and  for  latitude.  Based  on  geodetic  quantities  the 
equations  for  the  north  range  and  for  the  east  range  Rj-  in  the  local  tangent 
plane  are 
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%  ~ 

X  -  RgA$g 

(la) 

re  = 

Y  =  RgCos§gAX 

(1-1 

Rg  = 

geodetic  earth  radius 

A$g  = 

($gs  “  $gi )  radians 

AX  = 

(X8  -  X!  )  radians 

*g  = 

( $  gi  +  $  ga  )  /  2 

$gj ,  $ga  are  the  geodetic  latitudes  and  XT  ,  X8  the  longitudes  of  point  1  and  2 
respectively. 

From  standard  text  booics  (for  instance,  reference  7): 

Rg  =  Req(l  -  e2sin3$g)  2 

where 


Req  -  equatorial  radius  of  the  earth 

e  =  eccentricity  of  the  earth  ellipse 
Equation  (lb)  is  a  valid  approximation  and  can  be  used  for  quick  estimates;  (la) 
should  not  be  used  because  it  is  not  a  valid  expression  and  leads  to  significant 
errors  for  values  above  200  feet.  In  this  context  a  few  feet  are  considered  as  a 
tolerable  error.  Equation  (la)  is  mentioned  here  because  it  is  occasionally  used 
for  quick-look  estimates,  with  A$g  and  Rg  being  readily  available.  The  error  is 
shown  for  some  typical  values  of  Rjyj  in  the  section  on  "Numerical  Examples"; 
values  of  X(  =  R^)  computed  from  (la)  are  designated  as  "approximation  1". 

Selecting  a  certain  earth  model,  the  equatorial  radius  Req  and  the  polar  radius 
Rp0l  of  the  earth  are  given  and  can  be  used  to  compute  the  eccentricity  e  which  is 
defired  by 


=  1  -  (Rpol/Req>2  (2) 

A  second  approximation  for  X  uses  the  mean  geocentric  earth  radius  Rc  and 
geocentric  latitudes  $c: 

X  =  RCA$C  (3a) 


whe  re 


A?r. 


* 


rC3 


* 

t  c  \ 
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(3a)  is  a  valid  approximation  of  X  and  is  adequate  for  a  wider  range  compared  to 
X  computed  from  (la),  provided  the  two  points  are  on  the  same  meridian. 


(5a) 


The  geocentric  latitude  $c  can  be  computed  from 
$c  -  arctan[(Rpol/Req)3tan§g] 

Equation  (5a)  is  exact  and  is  found  in  reference  4.  $c  can  also  be  obtained  by 
applying  the  small  angle  approximation  to  tan($c  -  $ g)  and  then  expanding  in  a  power 
series.  This  is  done  in  appendix  A  and  the  first  two  terms  of  the  expansion  are 

$c  -  $g  -  -  0.  19390737sin2§g  -  0. 00131249sin2$gsin3  $g  (5b) 

where  $c  and  $g  are  in  degrees. 

Two  other  forms  of  this  series  are  given  in  references  2  and  3,  One  can 
derive  the  series  of  reference  3  by  manipulating  equation  (5b)  so  that  the  second 
term  in  (5b)  disappears  for  one  particular  value  of  $g.  This  is  described  in 
appendix  A.  Choosing  33“  for  this  particular  value  of  $g  yields: 

$c  -  $g  ^  -  0.  19429670s:n2$g  +  0.00038933(1-3. 371184sins$g)sm2$g  (5c) 

(Corresponding  to  equation  (A7)  in  appendix  A). 

If  one  neglects  the  second  term  of  (5c)  one  finds  an  improved  approximation  of 
($c  -  $g)  compared  to  a  one-term-only  solution  of  (5b).  This  improvement  is  ob¬ 
tained  of  course  for  the  specific  value  of  5g  and  to  some  degree  also  for  values  of 
$g  in  the  neighborhood  of  the  specific  value.  Since  the  mean  geodetic  latitude  at 
WSMR  was  33“  or  close  to  33“  for  those  target  ranges,  which  were  evaluated,  this 
number  was  considered  adequate  as  basis  lor  the  one-term  approximation. 

-  $g  -  -  0.  19429b70sin2$g  (~d) 

which  follows  from  (5c). 

The  values  of  X  and  Y  from  equations  (3a)  and  (3b)  can  now  be  computed  pro¬ 
vided  the  geode  to  latitudes  $g  are  given.  Examples  are  discussed  :r.  the  sect-on 

Numerual  Examples  .  Equation  (cd)  was  used  to  compute  ?  in  cquatior  (3a)  ard 
m  (3b)  and  the  approximations  for  X  and  Y  found  ir  this  way  are  referred  to  as 

approximation  2".  Th*  errors  caused  by  the  one-term  approximation  are  plotted 
,n  ticure  2  lor  >0*  and  3  3*. 

(n  pra.  tv  e.  another  procedure  >,s  more  frequently  used  to  estimate  the  length 
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of  an  arc  of  a  meridian  or  a  parallel,  and  this  estimate  in  turn  is  used  for  approxi¬ 
mating  north  or  east  components  in  the  local  tangent  plane.  This  is  the  method  of 
using  precalculated  scaled  factors  of  distance  in  feet  per  angular  degrees  or  per  arc 
minute.  These  scale  factors  can  be  derived  from  a  spheroidal  model  of  the  earth 
and  computed  as  functions  of  geodetic  latitude.  In  figure  3  and  tables  1  and  2,  these 
factors  are  presented  both  for  arcs  in  latitude  and  longitude  based  on  Clarke's 
spheroid  of  1866,  They  were  found  by  converting  a  one-arc -minute  difference  in 
latitude  or  longitude  respectively  to  tangent  plane  coordinates.  As  a  result,  the 
following  expressions  can  be  used  for  a  geodetic  latitude  of  33  to  estimate  the 
north/south  (X)  and  east/west  (Y)  coordinates  in  a  tangent  plane  system. 

X  =  6064. l£$g  (6a) 

Y  =  6093. 2AXcos33 ’  (6b) 

where  A$g  and  AX  are  in  minutes  of  arc,  or 

X  =  363848*  A$g  (6c) 

Y  =  365592AXcos33*  (6d) 

=  30661 1 AX 

where  A$g  and  AX  are  in  degrees. 

It  is  emphasized  here  that  equations  (3a),  (3b),  (6a),  (6b),  (6c)  and  (6d)  are  valid 
approximations  only  if  the  two  points  of  interest  are  located  on  the  same  meridian, 
or  or.  the  same  parallel,  respectively.  The  error  caused  by  not  being  on  the  same 
meridian  or  parallel  depends  on  the  value  of  the  angular  difference.  This  limitation 
is  related  to  the  subject  of  the  next  section, 
b.  Spherical  Ea  rth  Approximation 

A  bitter  approximation  is  obtained  bv  using  a  spherical  earth  model  with 
a  radius  equal  to  the  geocentric  radius  wh'ch  is  computed  as  the  mean  of  the 
gnu  cr.tr i  radw  of  the  two  particular  points  of  the  surface  of  the  earth.  The  ex¬ 
pressions  for  the  north  and  east  components  of  the  distance  between  two  points  in 
the  1  oval  tangent  plane  are  then  derived  from  figure  4  as  follows: 


8 


N 


! 


J 


Point  i 


(Rc+lfc  )cin*a  ^ 

:  /. 

(Rc+he  )co s$a  cosAX 


^*N 

(Rcth^  )c’o»#a,co8AXUn*1 


Rc+hi 


%  ertical  plane 
through  \x  -meridian 


AX 

•V 


(R_‘hp  I'  Oaij 


K. 

j 

■s 

Point  > 


r 


{  R  .  •  h-*  i  ‘0  5  l» 


1  Foirt  i 

i 

I 

i 

/ 


Equatorial  plar- 


Sarr  4.  Sph^r  .ial  Earth  Modrl 


l 


TABLE  ’ 


FEET  PER  DEGREE  AND  PER  ARC  MINUTE  AS  FUNCT'ON 
OF  LATITUDE  FOR  CLARKE'S  l«66  SPHERO’D 


Degree 

Feet  per 

Degree  of 
Geodetic  Latitude 

Fi  *  t  p«.  >* 

At  Mlr.ut, 
Geodet.t  Lat. 

0 

362, 760 

6046. 0 

10 

102, 874 

6  047.  6 

20 

363,  192 

6053.2 

30 

363, 684 

fcCt  1.4 

33 

363, 848 

6064.  i 

40 

36-,,  <190 

6r  ' ! ,  5 

50 

364, 632 

6082.2 

60 

365,  538 

6C92. 3 

70 

366,  036 

fc  1  Co.  o 

80 

? 66  .  3t 0 

‘  -06. ( 

Q0 

3c* , 4'4 

6  .C\  9 

Table  II 


FEET  PER  DEGREE  AND  PER  ARC  MINUTE  OF  LONGITUDE 
FOR  CLARKE'S  1S66  SPHEROID 

$ 

Degree 

Feet  per 

Degree  of  Longitude 

Feet  per 
Minute  of  Long: 

0 

365, 225 

6087. 1 

10 

359, 713 

5995.2 

20 

343,334 

5722.3 

30 

316, 562 

5276.0 

33 

306, 611 

5110.2 

40 

280, 170 

4669. 5 

50 

235, 230 

3920.5 

60 

183, 078 

3051.3 

70 

125,  289 

2088.  1 

80 

63,690 

1060.6 

90 

0 

0 
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It  can  be  seen  in  this  figure  that  Rj^  is  the  base  and 
(Rc  +  J^)sin$2  -  (Rc  +  hg  )cos$2  cosAXtar.^ 
is  the  hypotenuse  of  a  right  triangle  with  the.  angle  between  base  ar.d  hypoten¬ 
use.  Note  that  the  angles  and  $2  represent  geocentric  latitudes  and  R  is 
the  mean  geocentric  radius: 

Therefore 

X  -  RN  =  cos^j  [(Rc  +  hg  )sin$2  -  (Rc  +  hg  )tan$!  cos$2  cosAX] 

where 

AX  —  Xg  —  X  ^ 

or 

RN  =  (Rc  +  hg  Hcosfj  sin$8  -  sin$:  cos $a  cosAX )  (7a) 

Setting  AX  -  0,  under  the  assumption  that  the  two  points  are  on  the  same  meridian, 
one  obtains  as  can  be  expected 

Rn  =  (Rt-+ hg  )sin($g  -  $x  )  (7b) 

Only  for  small  angular  differences  ($2  -  )  and  1^  r  0  does  this  equation 

yield  approximately  the  same  results  as  equation  (3a)  which  indicates  the  limit¬ 
ation  of  the  circular-arc  approximation.  The  east  component  Rg  is 

Y  =  Re  =  (Rcg+  bfe  )cos$s  smAX  (7c) 

Where  Reg  is  the  geocentric  radius  of  point  2. 

Equation  (7c)  yields  an  approximation  for  Y  with  errors  depending  on  the  value  of  hg  . 
The  exact  value  of  Y  is  found  from 

Y  -  Re  =  (Rcacosfg  +  cos $g,s  )sinAX  (?d) 

It  is  noted  here  that  the  expression  for  the  north  component  (7a)  is  an  approxi¬ 
mation  because  of  the  use  of  a  spherical  earth  model  to  represent  an  elliptical 
earth;  but  (7d)  is  an  exact  expression  for  the  calculation  of  the  east  component  of 
the  range.  It  is  noteworthy  also  that  (7d)  uses  the  geocentric  radius  at  point  2. 

RCs  ,  instead  of  the  mean  geocentric  radius  Rc  which  is  used  m  (7a). 

In  many  practical  cases,  equation  (7a)  is  accurate  enough  and  its  simplicity 
makes  computation  possible  on  a  small  programmable  electronic  desk  calculator. 

An  example  of  the  program  which  was  prepared  for  the  Marchant  1016  PR  is  listed 
in  appendix  B.  The  program  consists  of  four  parts  each  of  which  is  recorded 
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}*c 

on  a  separate  magnetic  tape  .  These  tapes  are  read  into  the  core  memory 
sequentially  and  require  four  quantities  as  manual  input: 

fgj  =  Geodetic  latitude  (in  degrees)  of  point  1 

$gg  ~  Geodetic  latitude  (in  degrees)  of  point  2 

AX  =  Xa  -  Xx  (  in  degrees) 

1^  -  Height  of  point  2  in  feet  (multiplied  by  a  scale  factor  of  10"8 ). 

Typical  values  of  X  and  Y  were  calculated  from  equations  (7a)  and  (7c)  or  (7d) 

respectively  and  the  results  are  discussed  in  the  section  "Numerical  Examples". 
When  the  approximated  geocentric  latitudes  $c  from  equation  (5d)  are  used  to  com¬ 
pute  X  and  Y  the  approximations  are  referred  to  as  "approximation  3";  when  the 
exact  value  of  $c  from  equation  (5a)  is  inserted  into  equations  (7a),  (7c)  or  (7d) 
the  values  of  X  and  Y  are  referred  to  as  "approximation  4".  With  this  arrange¬ 
ment  the  approximations  are  easy  to  distinguish.  Those  designated  with  higher 
numbers  should  produce  more  accurate  results  with  respect  to  the  earth  spheroid. 

A  simple  relationship  (7e)  can  be  used  to  approximate  the  Z  component  of  the 
vector  connecting  point  1  and  2: 

Z  ~  [(Rea  +  ha)3  -  X3  -  YSP  -  (Rcl  +  1^ )  (7e) 

X  =  rn 

Y  *  re 

Note  that  (?e)  is  not  an  exact  equation  since  the  heights  hj  and  hg  do  not  form  a 
straight  line  with  Rcl  and  Rca  respectively;  however,  this  approximation  is  accu¬ 
rate  enough  for  the  applications  considered  here.  For  higher  accuracy  require¬ 
ments  the  Z  component  also  referred  to  as  R^»  should  be  derived  from  the 
spheroidal  earth  model. 

c.  Spheroidal  Earth  Model 

The  most  accurate  estimate  of  position  differences  is  obtained  by  applying 

a  spheroidal  or  ellipsoidal  model  of  the  earth  in  the  coordinate  transformation 

process.  Frequently  used  models  are  Clarke's  spheroid  of  1866  and  Hayford's 

*  With  minor  manipulations  explained  in  appendix  B  the  program  can  be 
recorded  on  3  tapes. 
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spheroid  of  1910,  also  called  the  international  spheroid.  Improved  earth  model 
parameters  were  recently  (1966)  derived  by  the  Smiths omar  Astrophvsic al 
Observatory  (SAO)  from  earth  satellite  data.  The  model  is  referred  to  later  i.r 
the  text  as  SAO  spheroid.  The  equatorial  and  polar  radius  for  these  three 
models  are  as  follows: 


—if 


Clarke's  spheroid  20,925*832  20,854,892 

Ha yford's  spheroid  20,926,470  20,856,010 

SAO  spheroid  20,925,738  20,855,576 

Derived  from  Rpol  and  R-eq  are  by  definition  the  ellipticity  E  -  1  -  Rpol/Req  and 

i 

the  eccentricity  e  =  [1  -  Rpol/Req)J  3  2,  which  are  listed  below  together  with 

e2 

and  the  frequently  used  term  j--- j 


_E _ e _ ef _ es  /  ( 1  -  e2  ) 

Clarke's  spheroid  1/295  .082271770  .0067686441  .0068147708 

Ha  yford's  spheroid  1/297  .08199218  .0067227183  .0067681701 

SAO  spheroid  1/298.25  .08182018  .0066945419  .0067396608 

Geodetic  survey  data  of  targets,  impact  locations  and  radar  sites,  etc.  ,  at  WSMR 
are  in  general  based  on  Clarke's  spheroid  of  1866  and  the  latter  is  vised  in  con¬ 
nection  with  numerical  examples  discussed  later.  There  are  two  ways  to  com¬ 
pute  the  north  and  east  components  of  the  range  between  two  points  from  giver, 
geodetic  data.  First,  the  "classical"  or  conventional  procedure  consisting  of 
the  conversion  from  geodetic  to  earth-centered  coordinates,  followed  by  a 
translation  and  rotation  of  the  coordinates.  Second,  relatively  simple  explicit 
formulae  can  be  derived  for  X,  Y  and  Z,  also  called  Rjsj,  Re  and  Rz  in  the  text. 


The  classical  method  is  described  as  follows: 

The  given  geodetic  coordinates  of  point  ore  and  two  are  converted  to  earth- 
centered  cartesian  system  (ecs)  which  is  a  left-handed  system.  The  coordinates 
of  point  two  are  translated  to  point  one  as  new  origin  which  yields: 
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xD  = 

^ecss 

"  ^ecsi 

YD  * 

Yecsa 

-  ^ecsi 

ZD  - 

^ecss 

^ecsl 

These  coordinates  are  rotated  three  times  to  obtain  the  final  set  of  cartesian 
coordinates  in  the  local  tangent  plane.  The  conversion  process  is  described  in 
more  detail  in  appendix  C.  Statements  of  a  FORTRAN  program  based  on  this 
procedure  are  included  in  appendix  D. 

The  other  procedure,  mentioned  above,  for  converting  geodetic  data  to 
local  tangent  plane  coordinates  is  described  next.  It  utilizes  in  a  straight¬ 
forward  manner  the  geometry  of  the  ellipse  representing  the  earth  and  of  the 
cartesian  coordinates  connecting  the  two  points  of  interest.  Similarly  as  for 
the  spherical  earth  model  in  figure  4,  point  2  and  straight  line  connections 
through  point  2  are  shown  in  figure  5  projected  into  the  vertical  plane  through 
the  meridian  on  which  point  1  is  located.  In  the  derivation  of  the  expression  for 
Rj^,  the  two  equations  of  the  vertical  and  horizontal  cartesian  coordinates  Ry 
and  Rpj  are  used  as  follows: 

Rv  =  Rg(l  -  c3  )sin$g  (9a) 

rH  -  Rgcos$g  (9b) 

Equation  (9b)  can  be  derived  from  figure.  5  readily  and  equation  (9a)  is  found  by 

inserting  (9b)  into  the  expression  for  tan$0 

Ry  b3 

tan$c  -  —  -  —  tan$„  (9c) 

RH  a  s 

which  yields 

RV  "  RgSin?g 

The  expressions  for  R^  and  R^  can  be  found  from  figure  5  which  shows  that  Rjq 
is  the  base  and 

H  =  (Ra  +  he  )sin$a  -c3(RBsin$8  -  R:  sin$j  ) 

-  (Re  +  he  )cos$a  cosAAtanSj  (9d, 
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is  the  hypotenuse  of  a  right  triangle  with  #j  being  the  angle  between  base  and 
nvpot«nuse.  Theref''»-e 

R-N  -  (Ra  +  he  )(cos#jSin$a  -  sin#j  cos$a  coaAX ) 

-  €2(R«  ain#8  -  Rj sinlj ) -osJj  (9e) 

As  can  be  expet  ted.  thi6  formula  is  similar  to  the  expression  for  Rj^  which  was 

derived  from  spherical  earth  (7c)  using  the  mean  geocentric  radius.  Since  (9e) 

is  based  on  the  elliptical  earth  model  utilizing  the  geodetic  radii  Rj  and  R®  and 

the  geodetic  latitude  #a  and  #8 ,  an  additional  term  must  be  applied.  This  term, 

the  second  part  of  (9e),  accounts  for  the  distance  between  the  centers  of  the 

radii  *3(Re  sinij  -  R.  sin#!  ).  The  expression  for  Rg  is  found  as 
*• 

Rz  =  Rj  +  1*  -  [{Rg  +  1^  )cos#a  cos^X /cos#!  +  Rpjtan#!  ]  (9f) 

The  equation  for  Rg  (9g)  is  equivalent  to  equation  (7c)  with  the  geodetic  radius 
Re  and  geodetic  latitude  #a  used  instead  of  the  corresponding  geocentric  data: 

Re  i  (Re  +  1*  )cos#a  sin^X  (9g) 

To  compare  the  two  procedures  described  above  for  the  computation  of  the 
LTP*c omponents ,  one  may  count  the  numerical  operations  (products,  divisions, 
dditions,  subtractions  and  squaring)  and  the  subroutine  entries.  Under  the 
assumption  that  In  either  case,  as  is  usual,  the  geodetic  coordinates  of  the  two 
points  ar<  giver.,  ;t  is  found  in  this  way  that  the  second  procedure  requires 

32  operations  including  7  trigonometric  functions 
ard  2  square  roots  to  compute  R^  ar.d  R^. 
and  ir  addition 

i  operations  to  i  ompute  R^. 

The  form  r  pro- 'dure  rt  quires  a  minimum  of 

41'  operations  r eluding  8  t* igonometr;*:  functions 
arc  2  square  roots  to  nomput.  and  Rg 

and  m  addition 

’  operations  to  compute  R^. 

Lo«a!  Tangent  Plane 
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A  double  precision  version  of  the  former  program  was  used  as  reference 
for  investigating  the  accuracy  of  the  two  programs.  It  was  found  that  the 
second  program,  the  new  version,  yields  slightly  more  accurate  results 
because  fewer  operations  are  involved.  The  difference  was  more  pronounced 
for  the  Z  component  than  for  the  X  component,  and  more  for  the  X  .ompom-nt 
than  for  the  Y  component.  Results  for  obtained  from  the  new  version  agree 
to  9  decimal  places  with  the  double  precision  results,  whereas  the  older 
procedure  leads  to  an  agreement  of  8  places. 

By  manipulating  the  matrix  product  in  equation  (C4)  of  appendix  C  repre¬ 
senting  the  classical  procedure  it  can  be  shown  that  (C4)  leads  to  equations 
identical  with  (9e),  (9f)  and  (9g). 
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SECTION  III 


NUMERICAL  EXAMPLES 


Accuracy  limitations  of  the  ciicular  arc  and  of  the  spherical  earth  approxima¬ 
tions  become  evident  from  numerical  examples  discussed  in  this  section.  The 
errors  AR^j  in  Rj^  are  presented  in  figures  6  through  10  which  are  made  by  the  4 
approximations  defined  as  follows: 


Approximation  1  circular  arcs  RgA$g 

Approximation  2  circular  arcs  RCA$C 

Approximation  3  spherical  earth  using  (5d)  for  ic 

Approximation  4  spherical  earth  using  (5a)  for  $c 

Those  values  of  Rjyj  served  as  a  reference  for  all  presented  errors  AR^  which  wer 
computed  from  equation  (9e)  based  on  Clarke's  spheroid  of  1866  for  various  ranges 
and  azimuth  angles.  These  reference  values  are  listed  in  tables  3  and  4  as  functions 
of  the  range  R^  (or  Rg)  and  of  the  azimuth. 

The  input  for  the  computation  is  assumed  to  be  given  in  geodetic  coordinates, 
latitude  $  and  longitude  X  in  degrees  and  height  h  in  feet  above  mean  sea  level.  For 

applications  considered  in  this  context,  it  is  adequate  to  carry  six  decimal  places 
for  $  and  X  and  to  use  a  rounded  integer  number  for  h  considering  a  one-foot 
accuracy  as  well  within  accuracy  requirements. 

Some  details  of  the  computation  which  apply  to  approximations  3  and  4  are 
clarified  first.  Neither  version  computes  the  mean  geocentric  radius  Rc  as  the 


mean  of  th«-  two  geocentric  radii  RC1  and  Rea  which  would  appear  necessary  from 
a  theoretical  standpoint;  instead  Rc  is  derived  from  the  mean  geocentric  latitude 
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Figure  9 


Approximation  3  and  4 


Azimuth  89*  to  91 v 


Figure  10 


AH  n  versus  Azimuth  Approximation  3  and  4 
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Rn  AND  Re  FOR  CLARKE'S  SPHEROID  .866  AS 
FUNCTION  OF  RANGE  AND  AZIMUTH 
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Tal:.t;  IV 

Rn  AND  R,e  FOR  CLARKE'S  SPHEROID  2  866  AS 
FUNCTION  OF  RANGE  AND  A  ZIMUTH 
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This  simplification  produces  ar  error  whi..h  is  insignifi -ant  c  ompared  to  other 
errors,  but  it  substantially  reduces  the  numerical  effort.  The  geo  ertrv.  earth 
radius  at  point  2,  Res,  is  used  to  compute  the  eas  +./'*■  r.:- 4  ;n*r»po-'’rl  Re  i'-nm 

Rg  =  Rea  cos$ca  sinAX  (l  Ob) 

This  expression  is  exact  and  values  of  Rg  computed  from  (10b)  a  rc  .dent:,  .a',  with 
those  derived  from  (9g)  if  the  exact  value  of  #<■*  is  used.  But  ever,  for  approximated 
geocentric  angles  $c 3  the  value  of  Rg  computed  from  (10b)  agreed  for  all  examples 
very  closely  (the  maximum  difference  found  was  3.4  feet)  with  the  value  of  R.g 
resulting  from  approximation  4.  Rg  values  are  therefore  omitted  from  further 
discussion. 

To  make  various  results  comparable,  the  assumption  was  made  in  all  cases 
that  the  center  of  the  range  between  the  aircraft  and  the  ground  target  was  at  3  V 
geodetic  latitude.  Basically,  two  parameters  have  been  varied  in  the  numerical 
examples.  These  are,  first,  the  length  of  the  range  vector  for  a  giver.  constant 
azimuth  and,  second,  the  azimuth  for  a  giver,  constant  range  vector. 

The  azimuth  Az  is  the  angle  between  north  and  the  range  vector,  counted 

positive  in  clockwise  direction.  It  is  computed  as  an  approximation  from  angu_ar 

data  (AX,  A$g)  instead  of  from  the  cartesian  coordinates. 

A.  =  arctan  . 

363848A  $g 

Six  different  azimuth  angles  are  considered  to  demonstrate  the  general  behavic*  of 
R^j  as  function  of  the  range  vector  length.  They  are  O',  45“  .  J  35'  89  9C  '  and  91  . 
For  these  angles,  Rj^  was  computed  from  equation  (9e)  as  tur.tr.or  o‘  the  *ar.ge: 
the  results  are  listed  ir:  Tables  3  and  4. 

As  can  be  expected,  R^  grows  nonlinear. .y  w  th  the  length  of  ore  of  the.  ve  tor 
coordinates  which  is  listed  in  nautical  mil*  s.  The  sign  of  is  positive  or  nega¬ 
tive  (same  as  the  sign  of  A?)  and  stays  so,  independent  of  the  vector  length.  An 
exception  to  this  rule  is  found  for  ar.  azimuth  slightly  larger  than  90°  as  car.  be 

seen  at  an  angle  of  91  .  (Table  4).  This  unusual  situatior  car.  be  explained  as 
follows: 
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The  sign  of  depends  in  this  spe  vial  ca  s*-  or  where  point  2  is  located 
with  respect  to  a  borderline  which  s  defi  .ed  by 

cosijsir.fe  cos§a  sir.fj  cosM  (11) 

This  equation  is  .our-d  from  (7a)  by  setting  R^j  0. 

If  this  equation  is  satisfied,  Rj.  vanishes;  :i  is  is  not  satisfied,  the  location  o‘ 
point  2  north  or  sc;th  determines  the  sign.  F  t  a  90  azimuth  therefore,  R^ 
is  a  relatively  small  positive  or  negative  quantity,  depending  or  whether  point 
2  is  west  or  east  of  point  1,  respe  t'.vely.  Increasing  the  vector  length  and, 
as  a  result,  also  increasing  Rg  for  the  case  ot  a  91°  azimuth,  point  2  changes 
its  position  from  ore  side  to  the  otht side  of  this  border  line  ar.d  R^  charges 
its  sign  accordingly.  This  situation  has  some,  effect  or:  ar.  application  dis¬ 
cussed  in  the  next  section. 

Considering  the  accuracy  of  the  cir  ;ula r  at,  approximation  first,  it  is  found 
from  figure  6  that  approximation  1  is  not  adequate,  for  ranges  RN  exceeding  0.  £ 
nautical  miles  and  approximation  2  i3  restricted  to  ranges  R^  up  to  33  nautical 
miles  provided,  ot  ourse,  that  the  two  points  are  located  on  the  same  meridian: 
for  larger  vaiu-.s  of  RN,  ar.  error  of  at  least  10  feet  is  made.  If  the  two  points 
are  not  on  the  sam-  meridian,  approximations  3  ard  4  can  be  used  over  a  com¬ 
paratively  vide  range  to  estimate  R^  wjth  tolerable  errors.  This  is  demonstrated 
ir.  figures  7  throegh  0. 

Inspecting  th-st  iigur*-s.  we  tird  that  as  a  gert-rai  trend  for  both  approxi¬ 
mations  of  RN  (approx  mat  ons  3  ard  4/  thf  -rror  £iR^  is  negative  and  its  absolute 
value  grows  as  a  rorl  .  *ar  .u'-'vr  with  th*  \ ».  -tor  length.  Its  value  varies  be  - 
tween  2C  ’eet  a*  d  1  feet  at  200  rauri  a?  m  Vs  der»-rd?r.g  on  th*  azimuth  argl* 
ar.d  reaches  mo;*  thar  tf'C  fe-t  «♦  :00  ru  al  mil*  s.  In  th«  neighborhood  of  90° , 
the  two  approximations.  ard  4,  y,«  d  almost  id*  rt- ...al  values:  the  differences 

are.  negligible  ior  pra  t.  al  purpose*,  i*tv>.;.  89  and  dl".  particular  interest 
is  th*  fa  t  that  for  the  azimuth  of  49'  th<  "poorer  approximation  (3),  repeegtrtir.g 

the  spher:  al  *  arth  model,  l*ads  to  a  valu*-  of  Rjyj  more  losely  resembling  the 


<0 


value  corresponding  to  the  spheroidal  earth.  This,  of  course,  is  an  effect  of 
the  bias  error  in  the  geocentric  latitudes  derived  from  the  trur.  at*  d  eqi.at.on 
mentioned  earlier. 

More  information  on  the  behavior  of  approximations  i  ard  4  is  fourd  >rom 
figure  10;  it  shows  for  a  constant  vector  length  of  about  150  nautical  miles  the 
errors  Rj^  as  functions  of  the  azimuth  angle.  The  maximum  rs  fourd  at  9C'\ 
From  this  diagram,  it  is  also  evident  that  within  some  inte rva„  between  10"  and 
90°,  approximation  3  yields  more  accurate  values  of  than  approximation  4. 
This  includes  the  case  of  45*  mentioned  earlier. 

It  is  noted  that  all  results  of  ARj^  shown  in  figures  fe  through  10  were  ob¬ 
tained  for  an  altitude  of  the  target  h,  equal  to  zero.  An  -yp*  rimer.t  showed  ’.har 
the  value  of  Ife  has  no  significant  effect  on  ARj^j.  Two  values  of  1^  were  con¬ 
sidered:  hj  -  0  and  ife  =  4000  feet  MSL,  which  is  approximately  the  tie  ration  of 
the  basin  of  White  Sands  Missile  Range.  The  effect  was  2  feet  or  less  vhen  Y.  , 
was  changed  from  0  to  4000  feet.  As  can  be  seen  in  equation  (9e),  1^  does  ha.e 
an  effect  on  Rj^  itself,  whereas  the  height  of  the  airplane  hj  does  not. 
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SECTION  IV 

IMPROVEMENTS  OF  THE  SPHERICAL  EARTH  APPROXIMATIONS 

Ar_  imp .rcv-.mer.t  of  th e  approximations  based  on  a  spherical  earth  model 

(app roximatior  s  3  and  4  in  the  previous  sections)  is  desirable  in  the  sense  that  Rj^, 
derived  from  a  slightly  modified  model,  approaches  more  closely*  the  talue  of  Rn 

which  is  based  on  a  spheroidal  earth  model.  This  is  desired  over  a  certain  in¬ 
terval  of  R^  values.  There  may  be  various  ways  of  achieving  this  improvement. 

As  it  has  teen  mentioned  earlier,  approximation.  3  may  yield  better  results  in 
this  respect  than  approximation  4;  the  improvement  is  caused  in  this  case  by  a 
bias  error  affecting  geocentric  latitudes.  Or  e  way  which  appears  practical  and 
efficient  or.  account  of  some  experiments  is  described  in  this  section. 

Considering  equation  (7a),  and  rewriting  It  as 

RN  =  Fx  Fa  (-  2a) 

vrher*- 

Fa  -  Rc  +  ha 

F3  cos}j  sin$3  -  ccs$3  sir.Jj  cosAX. 

it  is  evident  that  a  systematic  or  bias  type  t  rror  either  in  Fa  or  Fa  or  in  both  may 
affect  the  -  alee  ol  Rj^.  Introducing  a  bias  in  a  controlled  way  may  lead  to  the  im- 
p'or-  m*  t  o  the  approximation.  This  has  beer  done  b;,  adding  a  correction  term 


F j  us:. 

g  eh**  following  equations: 

F, 

(Rc  -  hg  .»  +  &R- 

AR 

ar.-0  Jb_  "I?*' 

Rn,  «c 

(12b) 

u  ?  ( RNr 

•  :*  O  gef.deC..  lat  e 

-  des  1  -  twe«  r  poirt  1  a-  d  2. 

s*  ie  ted  .  o:  vhi  .h 

1  V  adjustment. 

la  to  be  -educed  to  ?ero 

Rc 

m-.a r  geo..-.-  tri  :  earth  -ad 

lu  s . 

^rNVj 

th-.  •  r  •  i  r  to  be  ad  u?"ted  n~ 

compensated. 

rN’c 

r.h-  r  o-th  ccmpe:  -r*.  of  tre 

to  be  adiu^-t-  d. 

total  rarge  Sr  which  is 
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The  following  numerical  example  explains  the  procedure  lor  uncling  me  \aiue  ol 
ARC.  The  intention  is  to  make  ARjj  zero  at  a  range  Rj^  of  200  nauti-.al  miles 
for  an  azimuth  angle  of  135*.  The  corresponding,  uncompensated  error  ARN 
for  approximation  3  is  found  from  figure  8b  as 
ARj^j  =  -135  feet 

To  compensate  this  error,  Rc  must  be  reduced  by  a  proportionate  value 

Rc 

ARC0  -  -135  -  feet 

RNo 


With  Rc  =  20,904,910  feet  from  equation  (4b)  for  33*  and  with  Rjyjc  - 
1,183,256  feet  from  table  3  for  200  nautical  miles  (A2  =  135*),  we  obtain 


ARc 


-  2,  367  feet 
*  200  na 

6*(Rn) 


ARco 

For  values  of  R^  other  than  200  nautical  miles 


ARco 


Af0 


For  instance,  for  RN  =  100  nautical  miles,  using  the  value  of  A#  =  1.666"  f’-om 
table  3  and  6%=  3.32  for  200  nautical  miles,  we  find 


4Rc  .  .  2367  LUg 

=  -  1 1 87  feet 

In  the  same  way,  additional  correction  factors  car.  be  determined  between  R^  C 
and  Rpj  =  500  nautical  miles. 

For  later  reference  (12b)  is  rewritten  as 


AR-  --  C,A«(Rn)  ( 

where  Cj  is  a  constant  which  depends  on  the  azimuth  and  the  range  for  wh.'ch  aR^ 
is  compensated.  For  one  particular  vector  length,  the  error  AR^  is  >  omp  si’t-d 
and  for  values  above  and  below  it,  ARjsj  is  reduced  to  som-  r.xtert.  Nume  a', 
vab'^s  of  ARf  were  used  in  connection  with  the  Mar  hart  Cal  ulator  Prog'am 
described  in  appendix  B  to  compute  improved  app rox.mauo:  s  o:  ird  th« 

suits  are  referred  to  in  this  report  as  approximation  5  .  Results  obtained  v\ 
this  error  compensation  are  shown  in  figures  11,  12a  12b.  !l  and  14?.  For  C ' 

*  These  figures  also  contain  the  "ur  reduced  »rro"’  AR*y-  38  obtained  from 
approximation  4. 
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azimuth  (figure  11),  the  error  adjustment  is  of  p  tactical  value  only  at  large  values 
of  R«S|(R^>300  nautical  miles),  simply  t-.ca;?»:  approximation  4  leads  to  compar¬ 
atively  small  errors  ARjvj  at  lower  Rj^  values.  For  azimuth  angle  s  other  than  0*. 
a  significant  improvement  is  found  if  Ry  (or  Rj>;  exceeds  J  0C  nautical  miles. 

A  similar  reduction  of  ARiq  as  for  45*  and  1J5*  (figures  12a  and  i?t)  v/i-  cl  tai.zed 
for  67.5*  azimuth  also  by  applying  equation.  (13).  This  equation  can.be  used  in 
general  for  reducing  the  error  except  near  90*  azimuth.  For  90*  itself,  a  very 
simple  and  effective  error  correction  was  achieved  by  adding  a  constant  amount  of 
ARC  -  98925  feet 

to  the  mean  geocentric  earth  radius  R,~.  The  result  is  shown  in  figure  14.  It  is 
noted  here  that  the  value  of  this  constant  and  of  those  listed  below  d-p <■:/'.  on  the 
mean  geodetic  latitude  which  in  the  case  considered  here  was  always  33  ' . 

The  special  situation  for  azimuth  angles  slightly  bigger  or  smaller  than  90* 
requires  a  slightly  different  adjustment.  Considering  first  the  case  of  an  azi¬ 
muth  slightly  bigger  than  90*  (e.  g.  91*),  the  value  of  ARC  was  computed  from 

4R‘  ■  (!'” 

and  ARc  used  again  to  adjust  the  factor  F-  in  (J  2a}  as 
Fj  -  R«  +  hg  *  A  R^ 

The  resulting  reduced  error  ARjq  is  shown  in  figz'e  13.  dr.  a  small  region  which 
is  close  to  186  nautical  miles,  equation  (*4)  cannot  be  used  because  the  denom¬ 
inator  becomes  zero.  To  circumvent  this  difficulty,  another  set  of  coefficients 
Co  and  Cj  must  be  chosen  *. 

In  the  other  case,  if  the  azimuth  is  slightly  below  90*  (e.g.  89* ),  a  different 
equation  is  used  to  compute  AR^ •  A  second  degree  approximation,  was  found 
adequate  for  the  error  adjustment: 

ARC  -  CjAf  ♦  Cj  A#?  (15) 


*  A  better  solution  is  probably  to  use  two  sets  of  coefficients  for  equation  (13) 
corresponding  to  the  positive  and  negative  values  of  Ry. 


3? 


For  the  examples  shown  in  figures  11,  12a,  12b  and  13,  the  following  coefficients 

have  beer:  used  in  equations^!  3j  ^4}and  (i  5^  ! 

Azimuth 

(degrees) 


G 

C, 

90.030  fefct/d«.gree 

45 

Ci 

= 

399.24  f*;et/degroe 

135 

<n 

r 

713.08  feet/deglCe 

89 

.  Qi 

r 

1,219,  609  feet/degree 

Cj  -  5,  6?1 , 402  feet/degree3 

91  C0  -  0. 0943  degree 

C,  -  100,305.9  feet 

For  an  azimuth  of  67.  5"  the  coefficient  Cx  v/as  713.08  loot/degree. 
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SECTION  V 


CONCLUSIONS 

The  study  of  methods  ar.d  of  numerical  oxpr**ssior-s  for  the  oonv e  -  ■:  or,  o; 
geodetic  data  to  local  tangent  plane  coordinates  l*ads  to  the  following  cor.- 
elusions : 

1.  The  equations  for  computing  the  north  and  east  components  Rf*  ard  Re, 
based  or.  a  spherical  earth  model  can  be  mechanized  on  the  programmable 
Marchant  desk  calculator  1016  PR  and  stored  on  a  minimum  of  three  magnetic 
tapes  with  a  storage  capacity  of  100  bits  each.  The  program  uses  an  approxi¬ 
mation  for  computing  the  geocentric  latitude. 

2.  Depending  on  the  azimuth  of  the  local  tangent  plane  vector,  errors  of 
up  to  37  feet  in  R^  are  caused  by  this  program  for  a  100  nautical,  mile  dis¬ 
tance  in  north  (or  east)  direction  using  the  spheroidal  eanh  model  as  refer¬ 
ence.  This  error  grows  rapidly  with  increasing  length  of  the  vector.  The 
error  in  Rg  is  negligible. 

3.  The  error  in  Rjsj  as  mentioned  can  be  -reduced  by  using  precalculated 
correction  factors  which  depend  or.  the  azimuth  of  the  local  tangent ;  ’are 
vector.  Further  investigation  is  needed  to  determine  the.  optimum  correction 
factors  as  function  of  azimuth. 

4.  The  effect  of  the  height  of  the  target  point  (point  2)  in  the  interval  i-om 
0  to  4,000  feet  (MSL)  is  lound  insignificant  with  r-sperl  to  the  error  ARn- 
Discrepancies  in  the  results  did  rot  ex  re ed  2  f—.t.  for  all  ranges  Jr.v*- §tf gated. 

5.  A  simplified  equation  is  derived  to  d-t •-rrc.it*.  R R£  ar  d  R%  based  cr. 
the  spheroidal  earth  model:  the  rum-  rival  e.rort  -or  solving  it  is  "duc*-u  com¬ 
pared  to  the  effort  associated  with  the  star  da  "1  procedure. 

6.  Because  of  the  reduction  of  the  numerical  effort,  the  results  ar.d  sim¬ 
plifications  presented  are  applicable  fo'  quick  look  postflight  » -valuation  and  to 
some  extent  for  real  tim»-  data  processing. 
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APPENDIX  A 


COMPUTATION  OF  GEOCENTRIC  LATITUDE 


Based  on  equation  (5a)  in  the  basic  report,  a  power  series  is  derived  to 
approximate  )c  and  a  truncated  version  of  this  series  can  be  used  to  compute  { 
for  angles  in  the  neighborhood  of  a  given  value  of  #g.  To  find  the  power  senes, 
equations  (5a)  and  (2)  are  inserted  into 


tan(*g  -  wc)  = 


tan#  a  -  tan#c 


a  1  +  tan$ctan#g 

which  yields,  after  manipulation  of  trigonometric  relationships, 

r3  *in2#0 

tan(#g  -  #c)  =  1-  - - , 

2  1  -  J  sin  #g 

Assuming  that  #g  -  #c  is  a  small  quantity  and  applying  the  binominal  expansion, 
one  finds  the  difference  #g  -  fc  in  degrees  as 

#g  -  fc  “=  122  ii  sin2#g{l  +  ea  sina  *g  +  e*  sin*fg  +  . . . ) 

Solving  for  #t  and  carrying  only  one  term  yields 

fc  y  #g  -  122  (*1  +  «1  sin3#  )sin2#g  (A2) 

tt  2  2 

One  may  rewrite  (A2'  so  that  the  term  containing  sinafg  disappears  for  or.c 
particular  value  of  #g,  and  that  the  equation  for  #c  '•*  of  the  form 

#c  =  #g  -  Klsin2#g  +  K,(l  -  K,  sira  #g)sir.2#g  ...  (A3) 

To  determine  the  va’ue  of  the  coefficients  K*  ,  and  K*  ,  one  ma-  thm  proceed 
at  follows:  assuming  a  certain  special  value  for  fg  called  #g  9  is  c  hosen,  or.e 
substitutes  this  value  into  equation  (A2)  and  finds 

K,  -  112  («!-  ♦  il  sir.3  # „  ) 

TT  2  2 
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The  two  following  equations,  derived  from  (A2)  and  (A3),  are  then  compared  to 
determent  the  coefficients  and  1^ 

$c  ~  $g  -  'Kj  -  Kg  )sir.2$g  -  K3 111s  Ig*in2fg 

$  ?  «  -  l®o  r  «1  +  cl  8ina*  >  ,102*-. 

1  8  tt  ^  2  2  8  -  8 

To  satisfy  both  equations,  .t  is  found  that 
k_  .  180  e4  a. 

^  —  r  *8,  * 

K3  =  1 /sin2 *g4  b 

Inserting  the  numerical  value  for  smafg,  for  ♦g  -  30'  equation  (A2)  is  changed 
to  be 


*g~ 


!££  +  —  ^  sir.2*ff  +  JL_  (1  -  4sina  #(j)s:.n2# 

"  \  2  8  y  8  rr  8  8 


g 


In  a  similar  way  for  -  33' 

$.  ?  i  -  IM  (SL  a  il  0.  29663168^  s:r2*„ 

c  8  n  V  2  2  J  8 

+  —  si  0.  29663168  (1  •  3.  371 1 84s  ra #_) 

TT  2  •» 

Applying  a  value  for  f  corresponding  to  Clarke's  spheroid  1866,  ur.aliy  for 

$  -  30‘ 

*g.  * 


(A4, 


'AM 


(At) 


$,  r  '  ('0.  19423S49  -  0.0C032812  (l  -  4sm3  *g)}  s>.2ig 

and  tor  fg  s  31 

it  r  9  -  [0.  !9429fc70  -  C.  0C0  38VJ3  (i  -M  371  1 84s: M  $g)}  s.r  > (A' 

U'.r.g  *  quat  or  s  An  1  a'd  (A7»  th«  geo  er.t  r:  lat'tude  *  an  he  i  »!  ulai*  d  a*  ..-at*  1  v 
enough  »;>r  many  purpos*  s  .r.  the  *  •  lghborhcod  o:  30  or  33  g«  od-  tu  !a»  lude. 

Eq  ui!  ur  i.V  1  :s  !ou:  d  .r.  r*:«  r*-r,  r  3  w  th  r-»urn*  d-o:.‘  ;e  r  t  s.  Il  nr  \ .  or* 

ti  rm  s  us*  1  *  .  g.  tor  }g  #  'l 


y  -  j. 


-  0.  .  942  5'-4^  s  n2«. 


(AH  * 


ar.d  tor  s  ' 


\  t 


0.  1  9424b  70  i  r.2V 


g  .  ~'S 

a  stil!  relate*  ’v  small  «  rrur  •  made.  It  is  dtp*  ted  in  figure  2.  Equation  (A9>  :* 

us*  d  .r.  th*  Mar  har»  .fit-  PR  >*r»put«  r  program  d*  *  r.bed  :r.  apperdi*  B. 


(Adi 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  THE  MARCHANT  1 016  PR  CALCl'.’ATOR 

Purpose  ol'  the  program  was  to  lomput*  Irom  g_-  ►  r  geod’-t.  da'a  (inrp.t; o  , 
latitude  and  altitude)  the  tangent  plar.e  ?oord  rates  R^  and  Rjr  of  fort  B  w.th 
respect  to  the  origin  at  point  A  bas-d  or.  a  spher:  al  earth  model.  A  m.r  :mum 
of  memory  space  of  the  calculator  had  to  be  used  without  s .grift  art  loss 
of  accuracy.  The  program  corsists  of  4  sections  ea.  h  re  orded  or  a  mag¬ 
netic  tape,  read  into  the  calculator  irom  th*  IOTA  tap#  ur.t. 

Tape  I: 

This  part  of  the  program  computes  the  geo.entr:c  lat:tud’S  i  for 
point  l  and  point  2  using  the  approximation 
*c  =  *g  -  0.  1Q4296?  sir24g  ' 

whe  re 

}g  geodetic  latitude  ir.  degree# 

and  it  computes  the  mean  geocentric  iat  rude  f  . 

Input: 

Ig  j  ,  ig  a  ,  (in  d#  greesl  an-  load#  d  mar  uallv  .r  *o  tN  k*.vi<  4rd 

-»  giater 

Output: 

»  ,  in  W:  »_  ,  ir.  W2  t  ( i  .  .  f  s  ..  2  -  W  > 

All  three  angles  m  degree#. 

Cod* 

»  W2  m  x  x  w4-  m  _r  w 

1  M  N  X  Wb  R  ;  !  .  s  ■  ;  .  SW'ABO‘‘  HX  1' 

TRX  -1  .  Tv'  X  .  ?942«*- “  *W1  R  ■  S  M 

WI  R  BO  of  •  T  +  2  VI  M  K  T  M  K  BX  00 
'  This  expression  is  derived  -.e  app»'d  x  A. 

*  *  Marchant  symbols  art  used  for  the  tv  «  «•>  ,rt  :•>*•  t h«  ?  •  ga**\ «  s.gr  I-  the 
code  list  above.  N  stand  for  nega#  .  #  gr  .  \U  r  ha  t  .».  *  f  tor  r'  '  '  %rd 
n»  gat:\  r  sign. 


4' 


Controls  set  before  input  values  are  read  into  the  keyboard: 

10  digits  v,  clear  all  registers,  rour.d-off  switch  or.. 

Explanations: 

j  is  inserted  in  the  keyboard  before  the  "run  '  key  is  operated  and  the 
value  of  is  printed  immediately  before  the  actual  computation  starts.  In¬ 

sertions  of  $g  3  follows  when  the  program  comes  to  the  first  stop. 

Tape  2: 

This  part  of  the  program  computes  the  cosines  of  ^  ,a  *  *  an<* 

Input: 

}Cfl  in  Wl,  $Ci3  in  W2,  fc  in  W3  (all  in  degrees),  in  degrees  manually 
ir.to  the  keyboard  register  at  the  firct  stop. 

Output: 

cjsJc,  ,  in  Wl,  ios}CiJ  in  W2,  coslc  in  WJ,  cos&a  in  W4 
C  od  i, : 

-  X  .01745129252  X  W*.  M  W5  l  M  N  X 

W6  R  r  1  r  S  1  ♦  S  W5  A  BO  47  BX  21  T  R  + 

Wl  R  BO  7o  W2  R  BO  d4  W3  R  BO  92  T  W4  M 

K  T  ,\l  W2  BX  00  T  M  W3  BX  00  T  Ivi  K  #  BX  01 
Controls  set  h*  fore  running  section  2  of  the  program- 
’0  dig  ts.  round-off  sw.t  -h  on,  se'e„t  N'  . 

This  pa  rt  ‘imputes  tru  .man  geocentric  earth  radius  Rc  and  the  north 
c  oo  r  'i  •  :.at«  Kv  >>:  point  2  u  t  urgent  p’.are  w  ith  point  i  as  the  :<•>  girt  based  or  »n«‘ 

*  (JuAtcfo  t  B  •  at  u  ,  *  *  1 . 


1  digit  or*  s;.ir.  av.»*.d.»  ar  err  r  -if  .-  hr?,  •*  h  i  '  h  o.  i  * *  ii  .'be  » unction 
*:ri'  .*  >■::•>  -t.  :  ’■'th  i  digits  or.lv.  f.»r  ;  0.01  degree. 


4‘ 


whore 


fc8/(l-cw)]i 


0. 0067686441 


20,925,  832 
2.5348783 


’47. 7401 


Not*-  that  Rv  ,  A,  R«.q.  and  Rjq  art  ail  s  .a!*  d  dov.r  bv  a  factor  of  Th- 

formulation  of  (Bl)  reduces  the  memory  *pa •  *  requirement  to  a  mir.  mum. 


(Rf  t  hg  )(5:r.}  aoca}v  ,  -  s:r.  t  .  .  os#c  3  ,os£*  i 


In  j:  ut; 


cosie.i  ir*  Wl.  .oslr.  j  in  W2  cos#  in  W3,  «.  <js£*  ir  W4,  1^  r.to  k*.v 


board  at  the  first  stop. 
Printed  Output: 


Rk;  in  feet. 


Cooe: 


W  3  X  -14'.  "401  .  T/  -  2 .  r  •  -  8  r  r  ' 

T  I  A  ;;  f  A  W '  R  X  -1  v  T  •  X  W  .  R  .  M  T  X 
- !  -  T  /  X  W2  R  X  W4  R  -  -  W  '  R  ,  IX  W5  R  f  fcL 

Cm.tr  ils  set  before  ruT^rt  the  programs: 
it  d’jclts,  rour.d-oll  3w  :■  h  or. 

To  compute  the  fast  omponent  R;  •:?»•  the  a  1  ulat.o-  of  the  »*•  * 

. .»  repeated  ui:ft£  tap*  ’  ara  tape  Z  pti^-inu.  l  iter  tap.  4  »  r-  ad 
the  a!,  alator  to  -ompui* 


t  B  -  -  o  J  t-  3 


whe  re  as  l«lor«  • 


[e2/(l-e2)]i 

B  i=«l  +  1 

e* 

RE  :  (Rc,a  +  he)’'osr  a  sin&X 

Input: 

cos},  )3  in  W2,  cosAX  in  W4,  he  into  keyboard  at  the  first  stop. 
Printed  Output; 

Rg  In  feet 
Code: 

W2  R  X  --  -147.7401  +  T  /  +  2.  5348783  4  T 
W3  A  K  #  A  W4  R  X  =  -1  +  T  /  X  W3  X  W2  R  =  #  K 
Controls  for  running  the  program:  10  digits,  round-off  switch  on. 

The  codes  listed  above  are  not  necessarily  the  shortest  possible  codes. 


Remarks  on  Program  Operation  With  Three  Tapes  Instead  of  Four. 

With  very  little  extra  effort,  tape  4  can  be  eliminated  and  both  components, 
R^j  ar.d  R^.,  can  be  computed  using  only  3  tapes.  Tape  3  has  enough  unused  space 

to  arry  the  program  for  omputing  Rg.  After  Rjvj  has  been  computed  tapes  1  and 
2  ha  .  f  to  he  rerun  to  generate  r  os§^>2  unless  the  value  of  this  quantity  was  saved 
at  the  end  of  th»>  program  2  run  by  printing  the  contents  of  W2.  Before,  running 
tape  3  aga:.n  cos}  2  must  be  transferred  into  W3  (and  also  kept  in  W2).  The 
<  ode  oi  the  add, Mona!,  second  part  of  tape  3  for  computing  R.£  is  listed  below; 

W  4  R  X  --  - !  +  T/XW3RXW2R:|K 
The  stop  •  od<  K  at  the  end  of  the  first  part  of  the  tape  3  program  (see  page  47) 
must  be  replai  ed  by  the  !irst  letter  of  the  second  part,  which  is  W.  If  this 
oncratior  is  performed  with  only  3  tapes,  then  the  result  of  Rj£  printed  at  the  end 
of  the  first  path  and  the  result  of  R^  printed  at  the  end  of  the  second  path  are 
disrega  rd*  d. 
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APPENDIX  C 


TRANSLATING  AND  ROTATING  GEOCENTRIC 
COORDINATES  TO  A  LOCAL  TANGENT  PLANE 


This  appendix  derives  formulas  which  art.  used  ir-  a  star  dard  computer 
program  for  coordinate  transformation  under  the  option  for  orn  ertir.g  geodet: 
to  local  tangent  plane  coordinates.  The  pertinent  FORTRAN  statements  are. 
listed  in  appendix  D.  First,  the  exnressions  lor  earth  centere.d  (er-s)  coordi- 
nates  (equation  8) 

Xecs  --  RccosfccosX  (Cia) 

Yecs  -  RccosfcsinX  (Clb) 

ZeCs  “  Rcsinfc  ( C  * »  ) 

are  manipulated  by  inserting  the  equation  for  the  geocentric  earth  radius  R 
(equation  4b)n'  and  replacing  in  the  expression  for  Rc  the  term  e2  / ( 1  -e2)  by 
(a3  -b3  )b3 


where 

a  = 
b 

We  find 

X  ecs  - 
and  using  (da) 

Xecs  -  a3  (a”3  +  b3  tan3  §g) " ^cosX 
In  a  similar  way  we  find 

,2 ,  l 3 . 


equatorial  radius 
polar  radius 


a(l  +  tar3§c)  2cosX 
b 


i 


i 


v 

1  ecs 


a3  (ac  +  b3tan3 $g)"2sinX 


b2  (a2  +  b3tar3 


i 


"ecs  ~  w  xg >  "“"'g  1 

Note  that  in  the  equations  (CZ).  the  coordinates  are  functions  of  the  geo 

detic  latitude  which  is  more  readily  available  than  the  geocentric  latitude 
in  equations  (Cl ). 
v  equation  in  the  basic  report 


"Z  a) 


(C2h) 
’  C  2  ■.  ) 
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II  th<  ait  tern  h,  whi  h  was  assumed  as  r>  ro  so  J a ?  ,  has  a  fjo’tr  '  alue  we 
obtair 

X  I" a?  (a3  +  tan3  $„)"  2  +  h  rosl'oj  »  osA 

t-  c  S  '  g  g  J 

Yecs  -  I  a*  (a*  +  b^tan2  $g)“2  »  hos$gjs-r>. 

Ze(.s  -  b2(a2  +  b3tanc  in  )~  2tan$g  +  h  sin.S'g  (C3) 

The  following  steps  are  then  performed:  First,  the  e.  s  coordinates  are  computed 
for  point  1  using  (C3)  and  the  new  coordinates  ser\e  as  coordinates  of  the  origin 
of  the  new  system.  The  ecs  coordinates  of  point  2  are  computed  and  then  trans¬ 
lated  to  the  new  origin: 


XD 

x  PCS  ,3 

“  Xecs , x 

yd 

= 

Yecs,  3 

"  Yecs,  x 

7-d 

- 

Zees,  3 

“  ^ecs j ! 

The  coordinates  X^,  Y^,  Z^  are  rotated  in  three  steps.  First,  and  Y^  are 
rotated  by  an  amount  of  (+Xx )  degrees  around  the  vertical  Zp  axis  so  that  the  X 
axis  is  in  the  vertical  plane  through  the  meridian  of  point  1.  The  new  coordin¬ 
ates  are  called  X,  Y,  Z.  Second,  Y  and  Z  are-  rotated  by  18QJ  around  X. 

Finally  X  and  Z  are  rotated  in  north  direction  by  an  amount  of  (90  +  ) 

degrees  around  the  Y  axis  which  points  into  east  direction.  As  a  result,  X  points 
north  and  Z  upwards  along  the  local  vertical  at  point  1.  The  3  rotations  are 
expressed  m  the  following  equation 
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Th<  final  result  represents  a  left-hand*  d  system  with  X  r<  p)a.  <-fj  by  Y 

by  Rjr  and  Z  by 


Not*  :  The  component  Yq  in  west  direction  n  (C4)  has  a  regal  s  gi.  '<>  ohta.n 
positive  values  along  the  Yecs  axis  from  negative  longitudes,  in  th<  or  re¬ 
sponding  FORTRAN  program  in  appendix  Q  the  second  column  o',  the  '5  x  3  matr.x 
is  negative  instead  of  Yq: 


rn 

re 

-sinfcosX  -sin$sir.X 

-sinX  cosX 

cos# 

0 

’*d" 

yd 

RZ- 

cosfcosX  ros$sinX 

sir.l 

■« 

ZD 

^1 


APPENDIX  D 


FORTRAN  STATEMENTS  OF  CONVERSION  PROGRAMS 
1.  Cor  f-ntlona1  Program 

The  FORTRAN  statements  listed  he  r.  art  an  ex>  *  rpt  o'  a  program  for 
coordinate  transformation  which  was  c  oded  for  a  variety  of  applications  by 
D.  Dir- kins  or.  and  D.  Walter  at  AFMDC,  HA  Fit  in  1969,  Symbols  have  been 
changed  to  aid  the  understanding. 

Before  execution  of  the  program  at  the  computer  (CDC  3600/3800),  one 
card  is  read  to  ;rs-  rt  the  origin  of  the  LTP  system  with  -jooi  dinates: 

Geod'-hc  Latittude  $g,i  PHI  in  degrees 

Longitude.  X;  LA  Mi  ir.  degrees 

Altitude  h2  -  Hi  in  feet 

The  rext  (.ard  contains  the  geodetic  c  oordinates  of  the  point  2: 

Geodetic  Latitude  $g,a  -  PH2  in  degrees 
Longitude  X5  LAM2  ir.  degrees 

Altitude.  hg  H2  in  feet 

Input  data  are  based  or  Clarke's  spheroid  1866. 

Data  (DTR.  0.  '  '4--  3292619),  (R.EQ  2092^832), 

(RPC5L  208c"4 1. 92) 

State  m<  rt 
Numb,  i _ 

!  TYPE  REAL  LA  Ml  ,  LAM2 

2  PH  PH  DTR  $  LAM  1  LAM1  'DTR. 

3  PH2  PH2  DTR  $  LA  M2  LAM2  DTR 

4  SL2  ■-  SINF(  LAMc  '  $  C  L2  COSF(PH2) 

r>  SP2  SINFf  PH2 )  $  CP2  C(JSF(PH2)  $  TP2  SP2/CP2 


Statement 
Numbe  r 


6  A1  -  SQRTF({RP<2L  TP2)  2  +  REQ ;  2) 

7  R  -  REQ  '  *2/A1  t  H2  CP2 

8  XECS2  R*CL2 

9  YECS2  =  -RASL2 

10  ZECS2  =  RP0LJ‘*2*TP2/A1  +  H2  'SP2 

11  A21  -•  -SINF(LAMl)  $  A22  r  -CQ)SF(LAMl) 

12  A13  -  CG)SF(PHl)  $  A33  SINF(PH1 ) 

13  All  =  A22*A33  $  A12  =  -A21*A33 

14  A31  =  -A13*A22  $  A32  =  A13*A21 

15  A01  ■=  SQRTF(fRP0L*A33/Al  3)**2  +  REQ  ■  "2) 

16  RO  =  REQ**2/A01  +  H1*A1  3 

17  XECS1  -R0+A22 

18  YECS1  = +(R0*A21) 

19  ZECS1  ■=  RPOL**2*A33/(A13*A01)  +  HI  A33 

20  DX  =  XECS2  -  XECS1 

21  DY  =  YECS2  -  YECS1 

22  DZ  =  ZECS2  -  ZECS1 

23  RN  =  A?  UDX  +  A12  DY  +  A1  3  DZ 

24  RE  -  A2UDX  f  A22  DY 

25  RZ  -  A3UDX  4  A  DY  +  A3  5'  DZ 

New  Program: 

DATA  (DTR  =  0.  01  7453292619),  (REQ  =  20925b <2> 

(E2  «  0. 006768644065) 


r 


*  TYPE  REAL  LAM1,  LAM2 

PH i  -  PH1-DTR  $  PH2  =  PH2DTR 

*  DL  -  (LAMl  -  LAM2)  DTR 

4  SI  S  INF  (PHI )  $  CL  -  CQ>SF(PH1  )  ^  T!  -•  SI /Cl 
32  -  S1NF(PH2)  $  C2  C0SF(PH2) 

SDL  -  SINF(DL)  $  CDL  :  C0SF(CDL) 

7  R L  RLQ/SQRTF(1  -  E2  SI  -2) 

8  R2  REQ/SQRTF(  I  -  E2  •  S2  **2) 

9  RN  --  (R2  +  H2 )  *  (CDS2  -  S1*C2  CDL)  -  E2*(R2*S2  -  R1*SI)»C 

LC  RE  (R2  +  H2)'VC2'SDL 

•  l  RZ  R1  +  Hi  -  (R2  +  H2)V,C2~CDL/C1  -  RN*T1 
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